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NAVAL  AIR  DEVELOPMENT  CENTER 
AERO  ELECTRONIC  TECHNOLOOV  DEPARTMENT 
WARMINSTER,  PENNSYLVANIA  18974 


28  Dec  1972 


TECHNICAL  MEMORANDUM  NADC-20203 :OBL/PMM 


Sub! ’ Monthly  absolute  humidity  probabilities  for  selected  marine 
locations 

Ref:  (a)  AIRTASK  A360360C/001B/3F32343604 

(b)  Naval  Weather  Service  Command  Summary  of  Synoptic 
Meteorological  Observations 

(c)  NAVAIRDEVCEN  Tech  Memo  NADC-20203 :PMM,  Mathematical  Model  of 
PLIR  Performance  of  19  Oct  1972 

(d)  R.  D.  Hudson,  Jr.  Infrared  System  Engineering,  John  Wiley  and 
Sons,  Inc.  (1962) 


I . INTRODUCTION 

Reference  (a)  assigned  NAVAIRDEVCEN  the  task  of  performing 
operational  and  system  analyses  and  state-of-the-art  technology  surveys 
and  projections  as  a first  effort  in  the  development  of  FLIR  (forward 
looking  infrared)  equipments  which  would  be  affordable  in  large  quanti- 
ties and  optimized  for  the  missions  of  single-place  attack  aircraft. 

One  phase  of  this  task  is  the  calculation  of  the  PLIR  performance 
characteristics  (resolution,  sensitivity,  field  of  view  and  frame  rate) 
necessary  to  meet  operational  range  requirements  as  a function  of 
atmospheric  conditions,  aircraft  altitude  and  target  size,  orientation, 
and  effective  thermal  contrast.  In  these  calculations,  one  of  the 
crucial  factors  is  the  attenuation  by  the  intervening  atmosphere  of  the 
infrared  radiation  emitted  by  the  target  and  received  bv  the  sensor. 

Although  the  composition  of  the  atmosphere  is  quite  constant  with 
regard  to  most  of  its  constituents,  it  is  quite  variable  over  time  and 
place  with  respect  to  water  vapor,  which  is  a principal  absorber  of 
Infrared  radiation  over  the  spectral  range  of  interest  (3  to  14  micro- 
meters). Numerous  studies  of  the  transmission  of  infrared  radiation  as  a 
function  of  the  quantity  of  water  vapor  in  the  path  have  been  reported 
in  the  literature  and  several  methods  of  computation  have  been  published 
which  enable  one  to  calculate  atmospheric  spectral  transmission  as  a 
function  of  absolute  humidity  and  path  length.  Unfortunately,  there  is 
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a dearth  (or  perhaps  a complete  lack)  of  statistical  data  giving  probabil- 
ities of  occurrence  of  given  values  of  absolute  humidity  and  probabilities 
that  given  values  of  absolute  humidity  will  not  be  exceeded  in  certain 
times  and  places.  The  purpose  of  this  technical  memorandum  Is  to  provide 
marine  weather  statistics  in  a format  that  will  enable  one  to  determine 
the  absolute  humidity  for  which  a FLIR  must  be  designed  if  it  is  not  to  be 
prevented  by  atmospheric  water  vapor  from  achieving  its  operational  range 
requirements  to  any  preselected  degree  of  probability. 

II.  METHOD  OF  CALCULATION 

Reference  (b)  provides  monthlv  statistical  weather  data  for  a 
large  number  of  maritime  sites  throughout  the  world  based  largely  on 
reports  from  ships.  The  primary  period  of  record  for  these  data,  during 
which  the  most  recent  80%  of  the  observations  were  made,  is  typically  of 
the  order  of  20  years:  the  overall  periods  of  record  range  from  about  60 
to  110  years.  A representative  page  from  reference  (b) , which  gives 
data  for  the  Boston  area  in  May,  is  reproduced  here  as  figure  1.  Of 
particular  interest  in  this  figure  is  "table  13"  which  gives  the  percent 
frequency  that  the  relative  humidity  fell  within  each  of  eight  intervals 
of  relative  humidity  while  simultaneously  the  air  temperature  fell  within 
each  of  12  intervals  of  temperature  over  the  recording  period.  If  one 
takes  the  values  of  relative  humidity  and  temperature  at  the  center  of 
each  respective  interval  as  being  representative  of  the  whole  interval, 
one  can  then  compute  a representative  value  of  absolute  humidity  for 
each  of  the  96  two-dimensional  relative  humidity  - temperature  intervals. 
Reference  (c) , in  which  a convenient  formula  was  developed  for  performing 
such  calculations,  states  that  if  Hr  is  the  relative  humidity  expressed 
in  decimal  form,  T^  is  the  air  temperature  in  #C  at  sea  level,  and  e is 
the  base  of  the  Napierian  system  of  logarithms,  the  absolute  humidity  at 
sea  level  H0  is  given  to  within  +3.6%  for  temperatures  between  0°C  and 
35 #C  by 

0.06  T 

H0  “ 0-93  e A cm  of  precipitable  water  vapor  per  nmi. 


Alternatively, 


0.06  T.  , 

Hq  - 5.02  Hr  e A gm/mJ. 

These  calculations  were  carried  out  for  the  example  of  Boston  in  May  and 
the  results  are  summarized  in  table  I.  These  values  of  absolute 
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humidity,  in  conjunction  with  the  percent  freouencies  of  occurrence 
given  in  figure  1,  were  then  used  to  construct  the  histogram  of  figure  2. 
The  cumulative  percent  frequencies  of  occurrence  of  given  quantities  of 
atmospheric  water  vapor  concentration  were  then  obtained  by  integrating 
over  the  histogram  to  obtain  the  curve  shown  in  figure  3.  This  curve 
may  be  interpreted  as  representing  the  probability  that  the  absolute 
humidity  in  Boston  will  not  exceed  any  given  amount  during  the  month  of 


May. 


The  foregoing  approach  suffers  from  two  avoidable  weaknesses. 

The  first  of  these  is  that  the  absolute  humidity  calculated  from  the 
mean  temperature  and  mean  relative  humidity  associated  with  any  given 
relative  humidity  - temperature  interval  does  not  correspond  exactly  to 
the  mean  absolute  humidity  in  that  interval  because  absolute  humidity 
is  not  a linear  function  of  temperature.  The  second  weakness  of  this 
approach  is  that  it  gives  no  indication  of  the  magnitude  of  the  uncer- 
taintv  in  the  computed  results.  Because  of  the  way  the  data  are 
presented  in  reference  (b) , there  is  an  uncertaintv  in  temperature  of 
+2F#  and  an  uncertainty  in  relative  humiditv  of  (typically)  +5%.  The 
uncertainty  in  the  computed  values  of  absolute  humidity  is  larger  than 
either  of  these  taken  separately  because  there  is  a possibility  that 
the  highest  values  of  relative  humidity  (in  any  given  interval)  occurred 
at  the  same  times  as  the  highest  values  of  temperature  (in  the  same 
interval)  and  that  the  lowest  relative  humidities  and  temperatures 
(within  any  given  interval)  also  occurred  simultaneously.  The  foregoing 
can  be  Illustrated  by  an  example  taken  from  figure  1.  Consider  the  most 
probable  interval,  for  which  45 °F  <_  T^  < 49 °F  and  90%  <_  Hr  100%.  The 
highest  possible  absolute  humidity  that  could  fall  within  this  interval, 
which  occurs  if  T^  » 49#F  * 9.4*C  and  Hr  ■ 100%,  is  8.8  gm/m^.  On  the 
other  hand,  the  lowest  possible  absolute  humidity  that  could  fall  within 
this  interval,  which  occurs  if  T^  «*  45°F  • 7.2°C  and  Hr  ■ 90%,  is  7.0 
gm/m^.  The  absolute  humidity  for  this  interval  with  its  maximum 
uncertaintv  is  therefore  7.9  +0.9  gm/m^.  In  a similar  manner,  the 
absolute  humidities  corresponding  to  the  other  two  "comers"  of  the 
interval  were  also  calculated  and  all  four  values  are  given  in  figure  4. 

To  present  the  cumulative  percent  frequencies  of  occurrence  of 
the  various  concentrations  of  atmospheric  water  vapor  along  with  their 
maximum  uncertainties  it  was  decided  to  perform  two  absolute  humiditv 
computations  for  each  interval:  one  for  which  the  highest  temperature 
and  highest  relative  humidity  were  used  and  one  for  which  the  lowest 
temperature  and  lowest  relative  humidity  were  used.  Accordingly,  two 
curves  were  computed  for  each  selected  geographical  location  for  each 
of  the  12  months  of  the  vear;  one  of  these  curves  tends  to  overstate 
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the  absolute  humlditv  and  the  other  tends  to  understate  it.  In 
addition,  a third  curve,  representing  the  locus  of  mean  values  of 
these  extremes  of  absolute  humidity  was  also  plotted  for  each 
site-month. 

Nineteen  maritime  sites  in  the  northern  hemisphere  were  selected 
as  being  representative  of  the  types  of  climate  in  which  future  naval 
operations  might  be  expected  to  take  nlace.  Table  II  lists  and  describes 
these  sites  in  addition  to  providing  a two-letter  code  svstem  for 
identifving  by  site-month  each  of  the  graphs  to  follow.  Figure  5 
shows  the  geographical  distribution  of  the  sites  on  a small  scale  map 
of  the  world. 

Data  in  the  form  shown  in  figure  1 were  used  for  each  of  the  228 
site-months  for  which  computations  were  performed.  In  the  interest  of 
providing  greater  accuracy  over  a larger  temperature  range  than  that 
afforded  by  the  simple  equation  used  earlier  in  this  technical  memo- 
randum, values  of  absolute  humidity  corresponding  to  the  various  values 
of  temperature  and  relative  humidity  were  obtained  by  interpolating  and 
extrapolating  data  given  in  a table  in  reference  (d) . Cumulative  fre- 
quencies of  occurrence  of  the  various  values  of  absolute  humidity  were 
computed  by  means  of  a Hewlett-Packard  9100A  Calculator  and  plotted  as 
individual  points  by  a companion  9125A  Calculator  Plotter.  Smooth, 
eye-fitted  curves  were  drawn  through  each  set  of  points  and  then  curves 
representing  the  means  of  the  extremes  of  absolute  humidity  at  each 
cumulative  frequency  level  were  hand  drawn.  These  curves  are  included 
herein  on  the  228  graphs  labeled  AA  through  SL.  On  all  the  graphs  the 
cumulative  frequency  of  occurrence  is  plotted  as  the  ordinate,  ranging 
from  0 to  100%  while  the  absolute  humidity  expressed  in  grams  per  cubic 
meter  is  plotted  as  the  abscissa.  For  most  of  the  graphs  the  abscissa 
values  range  from  0 to  35  gm/m^;  however,  for  those  site-months  which 
very  high  absolute  humidities , this  scale  ranges  from  10  to  45 


Because  reference  (b)  gives  percent  frequencies  rounded  off  to 
the  nearest  0.1%,  the  maximum  values  of  the  cumulative  frequencies 
usually  do  not  equal  exactly  100%.  There  may  be  some  bias  in  the  original 
data  because  of  a preference  of  the  reporting  ships  to  avoid  extremes  of 
bad  weather  and  because  of  a preference  to  work  during  the  day  rather  than 
at  night. 

III.  CREDITS  AND  ACKNOWLEDGEMENTS 

The  task  of  locating  sources  of  data,  establishing  a method  of 
computation,  programing  the  Hewlett-Packard  calculator  and  performing 
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the  majority  of  the  computations  was  carried  out  by  Mr.  Herald  B.  Levin, 
a Temple  University  graduate  student,  who  was  a NA'TAIRDEVCEN  summer 
employee.  The  remaining  computations  were  performed  by  Miss  Nancy  E. 
MacMeekin  and  the  drawing  of  the  684  curves  was  the  work  of  Mr.  Glenn 
Jadney,  a student  trainee  from  Northeastern  University.  The  authors 
gratefully  acknowledge  the  contributions  of  Miss  MacMeekin  and  Mr.  Jadney. 
Special  appreciation  is  extended  to  personnel  of  the  Naval  Weather 
Service  Command  for  providing  the  relative  humidity  - temperature  data 
and  for  contributing  to  its  understanding  in  a number  of  lengthy  tele- 
phone conversations. 
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FI  CURE  2,  PERCENT  FREQUENCY  OF  OCCURRENCE  OF  GIVEN  VALUES 
OF  ABSOLUTE  HUMIDITY  FOR  BOSTON  IN  MAY 


Frequency 

of 

Occurrence 


Absolute  Humidity  (gm/m^) 


FIGURE  5.  GEOGRAPHICAL  DISTRIBUTION  OF  THE  SITES  FOR 
WHICH  ABSOLUTE  HUMIDITY  DATA  ARE  PRESENTED 


I"  VALUES  OF  ABSOLUTE  HUMIDITY  (gm/m  ) CORRESPONDING  TO 
VARIOUS  RELATIVE  HUMIDITY  - TEMPERATURE  INTERVALS 
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TABLE  II 

1 

1 IDENTIFICATION  OF 

SELECTED  SITES 

AND  MONTHS  BY  CODE  SYMBOLS 

FIRST  CODE 
LETTER 

LATITUDE 

LONGITUDE 

SITE  DESCRIPTION 

— 

A 

56  N 

151-157  W 

Kodiak,  Alaska  (NE  Pacific  Ocean) 

B 

45-47  N 

53-56  W 

Argentia,  Newfoundland  (Atlantic 

Ocean  SE  of  Newfoundland) 

C 

42  N-Coast 

66  W-Coast 

Roston,  Mass.  (Atlantic  Ocean  E of 
central  Mass.) 

D 

39.6  N 

129.0  E 

Wonsan,  North  Korea  (Sea  of  Japan) 

E 

36-38  N 

Coast-126  W 

San  Francisco,  Calif.  (Pacific 

Ocean  S and  W of  central  Calif.) 

F 

37.9  N 

25.1  E 

South  Aegean  Sea  (E  of  Athens , 

Greece) 

G 

36.0  N 

3.4  W 

Malaga,  Spain  (W  end  of  Mediterranean) 

H 

35.8  N 

133.6  E 

Matsue,  Japan  (S  end  of  Sea  of  Japan) 

1 

I 

32.2  N 

33.3  F. 

Port  Said,  Egypt  (SE  Mediterranean) 

J 

26.9  N 

50.3  E 

NW  Persian  Gulf  (between  Saudi 

Arabia  and  Iran) 

K 

25.1  N 

57.8  E 

N Gulf  of  Oman  (between  Iran  and 

Oman) 

L 

23.0  N 

67.7  E 

Karachi,  Pakistan  (N  Arabian  Sea) 

M 

20.4  N 

158.2  W 

Hawaiian  Leeward  (Pacific  Ocean  S 
of  Oahu) 

N 

18-20  N 

74-76  W 

Guantanamo,  Cuba  (Caribbean  S of  E 

Cuba) 

0 

17-22  N 

110  E 

S China  Sea  Area  VII  (E  of  North 

Vietnam) 

P 

k. 

17.8  N 

85.3  E 

Vlshakhapatnam,  India  (Bay  of 

Bengal  E of  S India) 
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TABLE  II  (cont'd) 


FIRST  CODE 

LETTER  LATITUDE  LONGITUDE SITE  DESCRIPTION 
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14.2  N 

73.0  E 

Panjlm,  Goa 
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(Arabian  Sea  W of  S 
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11-14  N 

111  E 

S China  Sea 
Vietnam) 

Area  I (E  of  South 
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7-11  N 

102  E 

S China  Sea 

Area  VI  (Gulf  of  Slam) 
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